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NOTE ON CERTAIN ITERATED AND MULTIPLE INTEGRALS 

By Wallie Abraham Hurwitz 

Disichlet* stated and used the following formula for inversion of the 
order of integration : 

/a fx pa ra 

dx I <f>(x, y)dy = / dy J <j>(x, y)dx. 

Volterrat makes use of a similar formula, which he calls "Dirichlet's 
principle." The truth of the theorem is obvious when the integrand is con- 
tinuous within and on the boundary of the field of integration. The object 
of the present paper is to justify the theorem when the integrand is allowed 
to become infinite in certain ways on the boundary of the field ; and to ex- 
tend it to space of n dimensions. Part II, which contains the generalisation 
to n dimensions, is entirely independent of Part I, and includes its results as 
a special case ; the method of treatment is, however, slightly different. 

While the results here stated, at least in so far as they relate to the case 
of two dimensions, may be deduced from a general theorem of de la Vallee 
Poussin,t the simple character of the reasoning here used seems to justify an 
independent treatment. 

I. 

A form of statement slightly different from that of Dirichlet and Volterra 
will be used, in which the integrals appear as. the generalisation of integrals 
mentioned by Schl6milch.§ 

Theorem. Letf(x, y) be continuous and limited in the region 

R: < x, < y, x + y < 1 ; 

* Crelle's Journal, vol. 17 (1887), p. 45. 

t Cf., for example, Annali di Matematica, vol. 25 (1897), p. 142. 

I Court d'analyse infinitesimal?,, vol. 2 (1906), §§ 73, 77. 

§ Analytische Studien, vol. 1 (1848), §19. Dirichlet, Llouville and other writers had also 
discussed special cases of the Integral here considered. Cf . the concluding paragraph of this 
paper. 
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and let < X, ft, v S 1 ; tfAew 

Jn ri — x 
f da; / x A -J y— >(1 - a; - yY~ 1 f(x, 2/)<ty 
o Jo 

= dy x*- 1 y-^l - a: - y)" _1 /(*, y)*K- (1) 

JO JO 

Let |/(x, y)\ ^ M. The integral 

/=/ ^/"'(l - x - y)— 1 f(x, y)dy 

is equivalent, by means of the change of variable y = (1 — x)s, to 

a^-^l - x)*+— ' f a"- 1 (1 - a)—i/[*, (1 - x)s2ds. 

The expression a:* -1 (1 — x)"+* ,_1 , which now stands outside the integral sign, 
is continuous for < a; < 1 ; the integral is readily seen to be uniformly con- 
vergent for ^ x ^ 1 , and hence to define a continuous function, whose ab- 
solute value is less than MB(fi, v). Thus the integral / represents a con- 
tinuous function for < x < 1 ; and 

\I\ ^ MB (ft, v) x*- 1 (1 - x)"+— \ 

Furthermore, I may be integrated between the limits and 1 ; for if we sub- 
stitute for it the above expression greater than \I\, we obtain the convergent 
integral 



MB (ft, v) f x*- 1 (1 - x)"+— 1 dx. 



Hence the iterated integral on the left of (1) is convergent. Since the two 
sides of (1) differ, essentially, only in having x and y interchanged, it follows 
that the iterated integral on the right of (1) also converges. It is evident that 
if we integrate over any intervals not reaching outside the intervals appearing 
in the integrals just considered, we shall still have convergent integrals. 

The field of integration is here a triangle bounded by the coordinate axes 
and the line x + y — 1. Draw lines inside this triangle, near its sides and 
parallel to them : 

x = 8, y = 8, x + y = 1 — 8. 
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Within and on the boundary of the smaller triangle thus formed the integrand 
is continuous ; integration may therefore be performed in either order with 
the same result: 

/l — 2« fl — x — t 

dx x*- 1 y*~ l (l - x - yy-i f(x, y)dy 

/— u ri—v—i 
dy / a*- 1 y-^l - * - yy-tfix, y)dx. (2) 

We may write 

/\— it n—x—t 
dx a^" 1 jr~ l (l -x- y)"- 1 ^, y)dy 

= ["dx r~V-iy"-i (1 - x - yy~ l f{x, y)dy - I x - I t , (3) 



where 



and 



/l —28 ft 

dx / s*- 1 jr- 1 (1 - x - y)"- l f(x, y) dy 

/— » n—x 
dx / a*- 1 y*- 1 (1 - * - y) v - l f(x, y)dy. 

Jl — X — t 

When y < 8, (1 - x — y)"~ x ^ (1 - x — S)— 1 ; hence 

/ ac*- 1 y-^l — as - y)»- x /(x, y)dy ^ itf**- 1 ^ - a; - S)"- 1 / y*~ l dy 
\Jo Jo 

= ^8*^-1 (i_a._sw ; 

/i— 28 | /-a 
dx / a5*- 1 y^ _1 (l —x — yy- l f(x, y)dy 
\ Jo 

/l — 28 
a*- 1 (1 - a: - S)"-i efcc 



*** 



M f 1 -* 
^—8* **-i(l -as- 8)— 1 dx 



(*) 
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Again, when y > 1 — x — h, y -1 ^ (1 — x — 8) 1 *— 1 ; hence 

Jri-x I 

I x»-'f->(l _ x -y)— i/(sc, y)<fy 
1— a— 8 I 

f (1 -x-y)'- l dy 

l—x—S 

M 

= ft &■ x*- 1 a - x - sy- 1 ; 

v v ' 

ri — M I ri—x 
therefore | 7 2 | — / dx / x A_1 y _1 (1 — x — y)"~ 1 f(x, y)dy 

M /* 1_M 

"■ x*- x (i -x-hy-^dx 



S — 8* / x*- 1 (1 - x - S)"- 1 dx 



= — B(X, v)fr(l -8)*+"-i. (5) 

From (4) and (5) we see that 

lim ij = 0, lim I 2 = ; 

and hence from (3) that the limit of the left-hand side of (2) is 

MC~* dx £~ x *~ 1 y " _1(1 ~ x ~ ^"~ l f( x > y"> dy 

= l dx l * xX ~ 1 y*~ 1 ( 1 - a:: -y) , "" 1 /0>y)<*y- (6) 

The two sides of (2) are obtained from each other by exchanging the roles of 
x and y ; it follows that also the limit of the right-hand side of (2) is 



M~ 



x*-*y»- 1 (l-x-yy-*f(x,y)dx. (7) 



If we now allow S to approach zero in (2), making use of the facts stated 
in (6) and (7), we obtain the theorem to be proved.* 

* It is in fact true that the double integral over the region B converges and is equal to 
the iterated integrals. The proof of this theorem, which would require separate treatment 
according as f(x, y) does or does not retain the same sign throughout M, is a special case of the 
more general considerations of Part II. 
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The form in which the theorem is used by Dirichlet, and by Volterra and 
other writers on integral equations, is obtainable from this result by means of 
the substitution as' = b — (6 — a)x. y' — a + (o — a)y. 
We find thus the 

Theorem. Let ifr(x, y) be continuous and limited within the triangle 
bounded by the lines x = y, x = b, y = a ; and let < X, p, v S 1 ; then 

I dx j(b - xy-^y - ay~ l (x - y)'~ l ^{x, y)dy 

Ja Ja 

= [ dy J (b - xy~\y - ay-\ x-y)— 1 ^ (x, y)dx. 

Ja Jy 

The preceding theorems may be stated for a wider class of functions. 
The equality of the results of iterated integration in the two orders in the 
interior triangle is sufficient for the validity of the proof as given. 

II. 

The preceding considerations are readily extended to the case of n varia- 
bles. We have first the following 

Theorem. Consider the n-dimensional region 



B 



[0<x, [»= 1, 2, ..-«]; 
(*! + Xj + • • • + x„ < 1 ; 



and the (w — 1) -dimensional region 

(0<x i [i=l, 2, ...n- 1]; 
B:\ 

(«! + X 2 + • • • + X B _! < 1. 

Let f(x u Xj, • • • x„) be continuous and limited in R, and let 

< \ u \j, • • • X„, X £ 1 ; 
then 



/x*." 1 ^- 1 



xjn-^l - xj - x 2 • - x^-^fdS 

1 — X,— *,— .. .— *. 

dS' 



? /"a^-ix*.- 1 • • • xj^-^l - x x - Xs- • • • - x„y-ifdx„. 
It should first be noticed that when f(x lt x 2 , • • • x„) = 1, the rc-fold 
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integral is known to be convergent, and its value is 

r(\ x )r(\,) . • • r(\„)r(\) * 
r(\j + \a + • • . + x„ + X) 

It will be convenient to write 

^-i j^-i . . . x\,-i(i — Xi — x 2 — • • • — a;,,)* -1 = 4>(x u ar 2 , • • • x„). 

Let |/| ^ ilf, and suppose at first that /SO. The integral 

1 — a. — — ar„ , 



'= / 4>dx„ 



subjected to the change of variable x„ = (1 — x 1 — • • • — x n _{)8, becomes 

a}*- 1 • • • ^ ,_1 (1 - *i • - *n-x) A+A - -1 

s*--^! - *) A-1 /[>i, • • • »»_!, (1 — *!—••• - x„_ 1 )s]ds. 



£ 



The integral which appears in the last expression converges uniformly through- 
out the region formed by R and its boundary, and therefore represents a con- 
tinuous function in this region ; the function denned by I is therefore of the 
same character in R as the given function <f> in R. f In R, I is continuous, 
and 

\I\ ^ MB{\, X„) a?'- 1 • • • afc.- 1 " 1 (1 - a, 3»- 1 ) A+A " _1 

Also we may integrate / over R, since the integral 

MB{\, X„) /si 1 " 1 - • • as>_y- x (l - as, a-O^" 1 ^ 

.As 

(which belongs to the special type just mentioned) is convergent. 

* Gonrsat, Cours cTanalyse, vol. I, §150. The proof there given under the hypothesis that 
the exponents are positive (so that the integral is proper) holds without alteration if each ex- 
ponent is greater than — 1. 

f If \ + X„ g 1, we take the function denned by the above integral with respect to a as 
the function corresponding to /, and X -f- X» as the number corresponding to X ; if X ■+■ ^n > 1 , 
we take the integral multiplied by (1 — %\ — • • • — x n _i) A +*»— l as the function corresponding 
to /, and 1 as the number corresponding to X. 
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Consider now the smaller fields of integration : 
(8 < Xi [i= 1, 2, • • • n] ; 
[x x + ajg + • • • + as» < 1 - S; 
_ (8<x ( [."=1, 2, ...n-1]; 
[x 1 + »j+ • • • +x„_ 1 < 1- 2S; 



10 < x 4 [i= 1, 2, • • • n 
«i + *j + • • • + s„_i < 1 



-1]; 
s. 

Since the integrand is continuous within and on the boundary of H s , 

1 — xi — ... — x n _j — S 

f <f>d8 = [d5 f <l>dx n . ( 8 ) 

As before, write 



, 1— x. — . . . — k , —t l—x, — ... — x„ 



<9) 



where 7 X = / dS I <l>dx n , 

J$ t Jo 

l_xl — . . . _ x 

and L 1= dS I <f>dx n . 

*« l-x, -...-«„_,-* 

When x»< 5,(1-^ «„)*-! S(l - x 1 as„-i - i)*" 1 ; 

hence 

J^dx n S itf^- 1 • • • x^~\l -x l x„_ x - $^- l Px x n - l dx n 



therefore 
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I, S ¥ 8K f x?- 1 • • • x^-\l -x, *„_,- «)*-i dS 

S *L BK fx^~ l ■ . • x*r, l_1 (l -SB, - ... - »._,- S)*-»tf£* 

Also, when 

x n > (1 -x a x„_i-S), x\. -^ (1 - *! as,-!-*)*-" 1 ; 

„1 — x, — ... — a., 

hence 



a — *i — . 

(l_ Xl *„)*-! 



.1 — X,— ... — x„ 

.., l dx„ 

l_x 1 -...-x„_ 1 -« 

iVf 
= _ jx^-i . . . ^'-'(1 _ x, *._,-*)».-* ; 

therefore 

7 2 S *f gA /" x*'" 1 • • • acjif-^l - x, *„_, - *)».-i dS 

S ^f 8> /x^'" 1 • • • x'lyi-^l - x t - • • • - as..! - S)*.-* d£ 

x r(x 1+ ■■■ + K-1 + K) K ' K ' 

From (10) and (11) we see that 

lim /j = 0, lim I 2 = ; 

6t=o s = o 

and therefore from (9) that the limit of the right-hand side of (8) is 

1 — X, — ... — *„_! 1 — X,— ... — *„_] 

l 1 "/ ^ ^ = / d ^/ * dXn " (12) 

* This integral is evaluated by means of the substitution x, = (1 — ff)Pfi[f = 1, 2, • • • n — 1] 
which reduces it to the special case already mentioned. 
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As for the left-hand side of (8), we note that since the integrand retains 
the same sign throughout R and since for a single set of regions R t whose limit 
is R, 

lim / $ds 
•-•/b, 

exists, therefore the integral over R converges and 



f 4>dS= lim J <t>d8. 



Hence we conclude from (8) and (12) that 



f 4>dS=fds[ 4>dx n , 
JR Jr Jo 



which was the theorem to be proved. 

Suppose now that,/ changes sign in R. It will always be possible to ex- 
press f as the difference of two functions which are nowhere negative, each 
of which is continuous and limited in R. * The theorem will hold for each of 
these functions, and hence for their difference /. 

We have thus reduced an n-fold integral to an (u — l)-fold integral of a 
simple integral. In the course of the proof it appeared that the function of 
n — 1 variables resulting from the simple integration has the same properties in 
R as the given function of n variables has in R. We may therefore repeat 
the process time after time, until we obtain an iterated integral in n variables. 
Furthermore, at each step the selection of the variable with respect to w'hich 
the simple integration is performed is evidently only a matter of notation ; we 
have therefore the 

Theorem. Let f{x l , x 2 , • • • x„) be continuous and limited in the region 

(0<x, [t= 1, 2, ... n]; 
R:\ 

and let <\ lt \$ • • • \„ £L 1; then the n-fold integral 

fx^- 1 ^- 1 . • • a^.-»(l - *! - x 2 - ■ • • - x n y~ l fdS 

* For instance, by writing/ = Jf — (M — /), where M is a constant greater than any 
value of / in JB. 
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converges, and, may be evaluated by iterated integration in any order; for ex- 
ample, 



Jit 



x^- 1 xfy- 1 • • • a5^» -1 (l — x 1 — x t — • • • — x n ) k - 1 fdS 



„i— *i 



= / x^~ x dx t / a;**- 1 dx 2 - ■ • I x^>- l (l — x l — x 2 — ---— x n ) K ~ l fdx n . 
Jo' Jo Jo 

(x'\ p i 
— \ ' and a slight change of notation, 

the theorem may be made to apply to the more general case in which the 
region of integration is 

< Xi \i= 1, 2, • • • n] ; 

I \aj \a.J \aj 

and the integrand is 

*-.*-i • • • <.-[! -(*)*-©- • • • -©'"]""><— "■*•>• 

where < X* S ft [z = 1, 2, • • • w] , < X S 1. 

Special cases of integrals of the types discussed in this paper were con- 
sidered by Dirichlet,* and shortly afterward by Liouvillef and Catalan. J 
The mosf general case (belonging to this type) mentioned by these writers is 
obtained by taking for f(x 1 , x 2 , • • • x n ) in this paper a function of x 1 + x 2 
+ • • ■ + x n . The problem discussed is always the evaluation of the integral, 
or at least its reduction to a simple integral, — all necessary considerations as 
to convergence and related questions being assumed. All these results are 
given by Schlomilch.§ References to later work of like character will be 
found in Encyclopadie der Mathematischen Wissenschaften, II A3, footnotes 
147-149. 

Harvard University, 

Cambridge, Mass., 

March, 1908. 

* Abhandlungen der Akademie der Wissenschaften zu Berlin, 1839, p. 61. 
t Journal de Mathematiques, vol. i (1839J, p. 225. 
X Journal de Mathematiques, vol. 4 (1839), p. 323. 
§ Loc. cit. 



